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The possibility of ferromagnetic and antiferromagnetic phase transitions in symmetric nuclear 
matter is analyzed within the framework of a Fermi liquid theory with the effective Gogny in- 
teraction. It is shown that at some critical density nuclear matter with the DIS effective force 
undergoes a phase transition to the antiferromagnetic spin state (opposite directions of neutron and 
proton spins). The self-consistent equations of spin polarized nuclear matter with the DIS force 
have no solutions corresponding to ferromagnetic spin ordering (the same direction of neutron and 
proton spins) and, hence, the ferromagnetic transition does not appear. The dependence of the 
antiferromagnetic spin polarization parameter as a function of density is found at zero temperature. 

PACS numbers: 21.65.+f; 75.25.+z; 71.10.Ay 



I. INTRODUCTION 



The spontaneous appearance of spin polarized states 
in nuclear matter is a topic of a great current interest 
due to its relevance in astrophysics. In particular, the 
effects of spin correlations in the medium strongly in- 
fluence the neutrino cross section and, hence, neutrino 
mean free path. Therefore, depending on whether nu- 
clear matter is spin polarized or not, drastically different 
scenarios of supernova explosion and cooling of neutron 
stars can be realized. Another aspect relates to pulsars, 
which are considered to be rapidly rotating neutron stars, 
surrounded by strong magnetic field. There is still no 
general consensus regarding the mechanism to generate 
such a strong magnetic field of a neutron star. One of 
the hypotheses is that a magnetic field can be produced 
by a spontaneous ordering of spins in the dense stellar 
core. 

The possibility of a phase transition of normal neutron 
and nuclear matter to the ferromagnetic spin state was 
studied by many authors Q-Qj predicting the ferromag- 
netic transition at p « (2-4) po for different parametriza- 
tions of Skyrme forces (go = 0.16 fm~^ is the nuclear 
matter saturation density). In particular, the stability of 
strongly asymmetric nuclear matter with respect to spin 
fluctuations was investigated in Ref. j2j, where it was 
shown that the system with localized protons can de- 
velop a spontaneous polarization, if the neutron-proton 
spin interaction exceeds some threshold value. This con- 
clusion was confirmed also by calculations within the rel- 
ativistic Dirac-Hartree-Fock approach to strongly asym- 
metric nuclear matter |10| |. Competition between ferro- 
magnetic (FM) and antiferromagnetic (AFM) spin order- 
ing in symmetric nuclear matter with the Skyrme effec- 
tive interaction was studied in Ref. , where it was clar- 
ified that the FM spin state is thermodynamically prefer- 



able to the AFM one for all relevant densities. However, 
strongly asymmetric nuclear matter with Skyrme forces 
undergoes a phase transition to a state with oppositely 
directed spins of neutrons and protons [l^ . 

For the models with realistic nucleon-nucleon (NN) in- 
teraction, the ferromagnetic phase transition seems to be 
suppressed up to densities well above ■'■'^ 
particular, no evidence of ferromagnetic instability has 
been found in recent studies of neutron matter [ig and 
asymmetric nuclear matter '17^ within the Brueckner- 
Hartree-Fock approximation with realistic Nijmegen II, 
Reid93, and Nijmegen NSC97e NN interactions. The 
same conclusion was obtained in Ref. where the 

magnetic susceptibility of neutron matter was calculated 
with the use of the Argonne two-body potential and 
Urbana IX three-body potential. 

Thus, the issue of appearance of spin polarized states 
in nuclear matter is a controversial one and models with 
effective Skyrme and realistic NN potentials predict dif- 
ferent results. From this point of view, it is interesting 
to attract another type of NN interaction and to com- 
pare the results for this NN potential with the previous 
results. Here we continue the study of spin polarizability 
of nuclear matter with the use of an effective NN interac- 
tion, namely, we utilize the effective Gogny force [l^[20l |. 
In addition, the reason to choose the Gogny interaction 
is as follows. It is known that the Skyrme interaction is 
a density dependent zero range NN potential. Its attrac- 
tive advantage is relative simplicity and successfulness 
in describing nuclei and their excited states. However, in 
many cases the finite range part of the nuclear interaction 
has the same importance as its density dependent zero 
range part |22| ■ This disadvantage of the Skyrme in- 
teraction is overcome by the Gogny interaction due to its 
finite range character. We will find the phase diagram of 
spin polarized nuclear matter with the Gogny interaction 
and will compare it with the results for the Skyrme and 
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realistic NN potentials. 

As a framework of consideration, we choose a Fermi 
liquid (FL) description of nuclear matter H^-jl^. We 
explore the possibility of FM and AFM phase transitions 
in nuclear matter, when the spins of protons and neutrons 
are aligned in the same direction or in the opposite di- 
rection, respectively. In contrast to the approach, based 
on the calculation of magnetic susceptibility, we obtain 
the self-consistent equations for the FM and AFM spin 
order parameters and solve them at zero temperature. 
This allows us not only to determine the critical density 
of instability with respect to spin fluctuations, but also to 
establish the density dependence of the order parameters 
and to clarify the question of thermodynamic stability of 
various phases. 

Note that we consider the thermodynamic properties 
of spin polarized states in nuclear matter up to the high 
density region relevant for astrophysics. Nevertheless, we 
take into account the nucleon degrees of freedom only, 
although other degrees of freedom, such as pions, hyper- 
ons, kaons, or quarks could be important at such high 
densities. 



II. BASIC EQUATIONS 

The normal states of nuclear matter are described by 
the normal distribution function of nucleons /kik2 = 
Tr ga'^^a^i , where k = (p, cr, r), p is the momentum, (t(t) 
is the projection of spin (isospin) on the third axis, and 
g is the density matrix of the system. The self-consistent 
matrix equation for determining the distribution function 
/ follows from the minimum condition of the thermody- 
namic potential [2^ and is 



/ = {exp(ro£ + n) + l} ' = {cxp(roO + 1} 



(1) 



Here the nucleon single particle energy e is defined 
through the energy functional of the system E{f) as 



.(/) 



dEjf) 



(2) 



In Eq. 1^, the quantities e and Y4 are matrices in the 
space of K variables, with 14^1^2 = YiTi^K.iK2 (''"1 = "-iP); 
Yq = l/T, Yin = -y^njT and Y^^ = -^ip/T being the 
Lagrange multipliers, /i„ and /ip being the chemical po- 
tentials of neutrons and protons, and T being the tem- 
perature. 

Taking into account the possibility of FM and AFM 
phase transitions in nuclear matter, the normal distribu- 
tion function / of nucleons and single particle energy e 
can be expanded in the Pauli matrices Ui and tu in spin 
and isospin spaces 

/(P) = /oo(p)'7oTo + /3o(p)cr3To (3) 
+ /03(P)(^0T3 + /33(P)(^3T3- 

£(p) = eoo(p)(ToTo + e30(p)CT3T0 (4) 
+ e03(p)croT3 + £33(p)CT3'r3. 



Expressions for the distribution functions /oo , /30 , /o3 1 /33 
in terms of the quantities e read 



(5) 



/o3 = -{n{uj+,+) ~ _) + n(tj_,+) - _)}, 

/33 = - n{uj+,^) - n(tj_,+) + _)}. 

Here n{uj) = {exp(ti;/T) + 1}^^ and 

^+,+ = Coo + C30 + "^03 + ^33, 
^+,- = Coo + C30 — Co3 — C33, 

^-.+ = Coo — Cao + Co3 — C33, 
t^-,- = Coo - Cso - Cos + C33, 

where 

Coo = ^00 ^ Moo J C30 = £301 

Co3 = £03 ^ Ai03, C33 = £33, 



(6) 



Moo 



M03 



The quantity being the exponent in the Fermi dis- 

tribution function n, plays the role of the quasiparticle 
spectrum. In the general case, the spectrum is fourfold 
split due to the spin and isospin dependence of the sin- 
gle particle energy £(p) in Eq. Q). The branches lo±.+ 
correspond to neutrons with spin up and spin down, and 
the branches lo±_^ correspond to protons with spin up 
and spin down. 

The distribution functions / should satisfy the norma- 
lization conditions 



/o3(p) = Qn - Qp^ ag, 
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(7) 
(8) 

^Xl-^3o(p) = £•! - ^?i = ^PTT' (9) 
/33(p) = (ftiT + evi) - (ftU + 6'pt) = ^^Tl- (10) 



Here a is the isospin asymmetry parameter, gi„|, g^i and 
Qv\ I 6vi neutron and proton number densities 

with spin up and spin down, respectively; g^ = p„| -|- gp-^ 
and gi — -|- are the nucleon densities with spin 
up and spin down. The quantities ^g\\ and Apt 1 play 
the roles of FM and AFM spin order parameters ■ 

In order to characterize spin ordering in the neutron 
and proton subsystems, it is convenient to introduce neu- 
tron and proton spin polarization parameters 



n„ 



Qn 



H„ 



Qp 



(11) 
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The self-consistent equations for the components of 
the single particle energy have the form [ill [l3 | 



Coo(p) = eo(p) + £oo(p) - A^oo, 6o(p) = £3o(p), (12) 

C03(P) = £03 (P) - C33(P) = £33 (P)- 



t-2 2 

Here eo(p) = -5;^ is the free single particle spectrum, 
mo is the bare mass of a nucleon, and eoo, £30, £03, £33 are 
the FL corrections to the free single particle spectrum, 
related to the normal FL amplitudes Uo(k), J73(k) by 
formulas 



£oo(p) = ^Et^o(k)/oo(q), (13) 
q 

£30(p) = ^E^lW/30('l)' 

q 

£03(p) = ^E^2(k)/03(q), 

q 

£33(p) = ^E^3('')/33(q)- 

q 



Further we do not take into account the effective ten- 
sor forces, which lead to coupling of the momentum and 
spin degrees of freedom [23, l27l l2^. and, correspond- 
ingly, to anisotropy in the momentum dependence of FL 
amplitudes with respect to the spin polarization axis. 

To obtain numerical results, we use the effective Gogny 
interaction. The amplitude of the NN interaction in this 
case reads 



2 

^(P,q) - toil + XoPa)g^ + ^M'(W^^ + B,P„ 

1=1 

(14) 

- H,Pr - M,P„A)e-(P-'^)''^'/^ 

where Pq- and _FV are the spin and isospin ex- 
change operators, and tg, xq, fii,Wi, Bi, Hi and Mi are 
some phenomenological constants, characterizing a given 
parametrization of the Gogny forces. In numerical cal- 
culations we shall utilize the DIS potential Using 
the same procedure as in Ref. it is possible to find 
expressions for the normal FL amplitudes in terms of 



Gogny force parameters 

2 

[/o(k) = eiopT + 27r3/2 J2 M?(2S, - 2H, - M, + AWi) 
1=1 

(15) 

2 

- 27r3/2 e-''"^'"^i|{2B, - 2H, - 4M, + Wi), 
1=1 

2 

Ui{k) = -2tog^{l - 2xo) + 2^3/2 ^ ^3^2B, - M,) 

i=l 

i=l 

2 

C/2(k) = -2toe^(l + 2xo) ~ 2^3/2 ^ ^3(2ij^ + M{) 

i=l 

-27r3/2^e-''^A<?^3(2s^_^_^^)^ 

i=l 

f/3(k) = -2tog^ 

i=l i=l 

Thus, with account of expressions Q for the distribu- 
tion functions /, we obtain the self-consistent equations 
(|12|l . H13() for the components of the single particle en- 
ergy 60 (p) , 60 (p) , 63 (p) , 63 (p) , which should be solved 
jointly with the normalization conditions (f7|)- (|10|l . deter- 
mining the chemical potentials Hoo^fJ-os, FM and AFM 
spin order parameters Ap||, Agi^^. Since the FL ampli- 
tudes in Eqs. p5|) contain two Gaussian terms, the self- 
consistent equations represent, in fact, the set of coupled 
integral equations, which can be solved iteratively using 
the Gaussian mesh points in the momentum space. 



III. PHASE TRANSITIONS IN SYMMETRIC 
NUCLEAR MATTER 

Early research on spin polarizability of nuclear matter 
was based on the calculation of magnetic susceptibility 
and finding its pole structure 4, 5], determining the onset 
of instability with respect to spin fiuctuations. Here we 
shall solve directly the self-consistent equations for the 
FM and AFM spin order parameters at zero temperature. 
In this study we consider the case of symmetric nuclear 
matter (p„ — Qp). 

The FM spin ordering corresponds to the case 7^ 
0, ^30 (p) ^ 0, Ag^i = 0, C33(p) = 0, and there are two un- 
known parameters /zqo, Ag-^^ and two unknown functions 
Coo(p),6o(p) (mo3 = 0,eo3(p) = as a consequence of 
isospin symmetry) . The AFM spin ordering corresponds 
to the case Ag^l ^ 0,63(p) / 0, Ag|| = 0,6o(p) = 
and we should find two unknown parameters noo^Ag^i 
and two unknown functions Coo(p)j C33(p)- 
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In the FM spin state of symmetric nuclear matter we 
have = Qp-\,Qni = Qpi and nucleons with spin up 
fill the Fermi surface of radius k2 and nucleons with spin 
down occupy the Fermi surface of radius fci , which satisfy 
the relationships 

^{kl~kl)^Ag^T, ^(kl + kl) = g. (16) 

Since at zero temperature there are no spin up nucleons 
with k > k2 and there are no spin down nucleons with k > 
ki, then, as follows from Eq. w+^+(fc2) = (jj+^_(fc2) = 
0, w_,+ (fci) :=0. 

In the AFM spin state Qni = Qpi, Qni — Qp^ and neu- 
trons with spin up and protons with spin down fill the 
Fermi surface of radius fc2 and neutrons with spin down 
and protons with spin up occupy the Fermi surface of 
radius fci, satisfying the equations 

^(fc|-fc3) = A^,u, ^{kl + kl)^e. (17) 

At zero temperature there are no spin up neutrons 
and spin down protons with k > k2 and there are no 
spin down neutrons and spin up protons with k > ki. 
Hence, as follows from Eq. ©, oj+,+ {k2) = a;___(A;2) = 0, 
^+,-(^i) = w_,+ (fci) = 0. 

In the totally FM polarized state we have Agi|| — 
Q, k2 = kp = (Stt^p)^/'^ and the degrees of freedom, 
corresponding to nucleons with spin down are frozen 
{ki = 0). For totally AFM polarized nuclear matter we 
have Aqii = g, k2 = kp, where k2 is given by the same 
expression as in the FM case, since now the degrees of 
freedom, corresponding to neutrons with spin down and 
protons with spin up, are frozen (ki — 0). 

Now we present the results of the numerical solution 
of the self-consistent equations with the DIS Gogny ef- 
fective force. The main qualitative feature is that for the 
DIS force there are solutions corresponding to AFM spin 
ordering and there are no solutions corresponding to FM 
spin ordering. The reason is that the sign of the multi- 
plier to{2xo — 1) in the density dependent term of the FL 
amplitude Ui, determining spin-spin correlations, is pos- 
itive, and, hence, the corresponding term increases with 
increase of nuclear matter density, preventing instability 
with respect to spin fluctuations. Contrarily, the density 
dependent term —2tog'^ in the FL amplitude U3, describ- 
ing spin-isospin correlations, is negative, leading to spin 
instability with oppositely directed spins of neutrons and 
protons at higher densities. Here the situation is similar 
to that with the Skyrme effective forces SLy4 and SLy5 
in strongly asymmetric nuclear matter [T^ . where anal- 
ogous behavior of the FL amplitudes Ui and U3 in the 
high density domain prohibits the formation of the state 
with the same direction of neutron and proton spins and 
leads to the appearance of the state with the oppositely 
directed spins of neutrons and protons at high densities. 
However, the results with the Gogny effective interaction 
are in contrast with the results of microscopic calcula- 
tions with a realistic NN interaction jT^ , predicting that 
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FIG. 1: AFM spin polarization parameter as a function of 
density at zero temperature for the DIS Gogny force and the 
SkM*, SGII Skyrme forces. 

the FM spin state is always preferable over the AFM one 
for all relevant densities, but is less favorable compared 
to the normal state. 

In Fig. n] it is shown the dependence of the AFM spin 
polarization parameter Ag^i/g as a function of density 
at zero temperature. The AFM spin order parameter 
arises at density g w 3.8^0 for the DIS potential. A 
totally antiferromagnetically polarized state {Ag^i/g = 
1) is formed at g « 4.3po- The neutron and proton spin 
polarization parameters for the AFM spin ordered state 
are opposite in sign and equal to 

n --H 

iin — lip — 

g 

For comparison, we plot in Fig. ^ the density depen- 
dence of the AFM spin polarization parameter for the 
Skyrme effective forces SkM* and SGII 0. It is seen 
that the results with the DIS potential are close to those 
with the SkM* potential (for the DIS force the AFM 
spin polarization parameter is saturated within a nar- 
rower density domain than for the SkM* force). 

To check the thermodynamic stability of the spin or- 
dered state with oppositely directed spins of neutrons 
and protons, it is necessary to compare the free ener- 
gies of this state and the normal state. In Fig. |21 the 
difference between the total energies per nucleon of the 
spin ordered and normal states is shown as a function of 
density at zero temperature. One can see that nuclear 
matter in the model with the DIS effective force under- 
goes at some critical density a phase transition to the 
AFM spin state. 

In Fig. 13 the difference between the total energies per 
nucleon of the spin polarized and normal states is decom- 
posed into two parts, the kinetic and correlation ones. In 
spite of increasing the kinetic energy per nucleon in the 
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FIG. 2: Total energy per nucleon, measured from its value 
in the normal state, for the AFM spin state as a function of 
density at zero temperature for the DIS Gogny force and the 
SkM*, SGII Skyrme forces. 
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FIG. 3: Total energy per nucleon, measured from its value 
in the normal state, for the AFM spin state as a function of 
density at zero temperature for the DIS Gogny force, decom- 
posed into kinetic and correlation parts. 

AFM spin state, the AFM spin state becomes thermo- 
dynamically preferable over the normal state due to the 
energy gain caused by medium correlations, mainly by 
spin-isospin correlations, leading to AFM instability of 
the ground state. 



IV. DISCUSSION AND CONCLUSIONS 

Spin instability is a common feature associated with a 
large class of Skyrme models, but is not realized in more 



microscopic calculations. In this respect, it is interesting 
to study the possibility of appearance of spin polarized 
states in nuclear matter, utilizing another type of NN 
potential. Here we used the finite range Gogny effective 
interaction, being successful in describing nuclei and their 
excited states. The force parameters are determined em- 
pirically by calculating the ground state in the Hartree- 
Fock approximation and by fitting the observed ground 
state properties of nuclei and nuclear matter. The anal- 
ysis based on the Gogny interaction shows that the self- 
consistent equations of symmetric nuclear matter have 
solutions corresponding to AFM spin ordering and have 
no solutions at all corresponding to FM spin ordering. 
This result is in contrast to calculations with the Skyrme 
interaction, predicting a FM phase transition in symmet- 
ric nuclear matter [ij. In the last case, the SkM*, SGII 
parametrizations were used, developed to fit the prop- 
erties of nucleon systems with small isospin asymmetry. 
However, if the spin ordered state with oppositely di- 
rected spins of neutrons and protons for the Gogny in- 
teraction survives in the domain of high isospin asymme- 
try, then this result can coincide with the calculations for 
the Skyrme interaction, predicting the appearance of this 
state in strongly asymmetric nuclear matter '12]. In the 
last case, the calculations were done with the SLy4, SLy5 
parametrizations, developed to reproduce the properties 
of nuclear matter with high isospin asymmetry. 

In a microscopic approach, one starts with the bare 
interaction and obtains an effective particle-hole interac- 
tion by solving iteratively the Bethe-Goldstone equation. 
In contrast to the Skyrme and Gogny models, calcula- 
tions with realistic NN potentials predict more repulsive 
total energy per particle for a polarized state comparing 
to the nonpolarized one for all relevant densities, and, 
hence, give no indication of a phase transition to a spin 
ordered state. 

It must be emphasized that the interaction in the spin 
and isospin dependent channels is a crucial ingredient in 
calculating spin properties of nuclear matter and differ- 
ent behavior at high densities of the interaction ampli- 
tudes, describing spin-spin and spin-isospin correlations, 
lies behind this divergence in calculations with effective 
and realistic potentials, from one side, and calculations 
with different types of effective forces, from the other 
side. Since our calculations with the Gogny interaction 
predict the AFM spin state as the ground state of sym- 
metric nuclear matter, this emphasizes the role of spin- 
isospin correlations in the high density region. Due to the 
antiferromagnetic spin polarization, some neutrons and 
protons with opposite spins, e.g., spin up neutrons and 
spin down protons, fill the Fermi surface with the larger 
radius and others, spin down neutrons and spin up pro- 
tons, occupy the Fermi surface with the smaller radius. 
When density increases, some neutrons and protons un- 
dergo spin-flip transitions from the inner Fermi surface 
to the outer one due to increase of spin-isospin correla- 
tions. The usual way to constrain the interaction param- 
eters of spin-dependent amplitudes is based on the data 
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on isoscalar [23 and isovector (giant Gamow- Teller) [S^l 
spin-flip resonances. However, it is necessary to note 
that in order to get robust results for the spin polariza- 
tion phenomena, these constraints should be obtained for 
the high density region of nuclear matter. Probably, such 
constraints can be obtained from the data on the time 
decay of the magnetic field of isolated neutron stars [sJl ■ 
In summary, we have considered the possibility of 
phase transitions into spin ordered states of symmetric 
nuclear matter within the Fermi liquid formalism, where 
the NN interaction is described by the DIS Gogny ef- 
fective force. In contrast to the previous considerations, 
where the possibility of formation of FM spin polarized 
states was studied on the basis of calculation of the mag- 



netic susceptibility, we obtain self-consistent equations 
for the FM and AFM spin order parameters and solve 
them at zero temperature. It has been shown in the 
model with the DIS effective force that symmetric nu- 
clear matter undergoes a phase transition to the spin po- 
larized state with oppositely directed spins of neutrons 
and protons, while the state with the same direction of 
the neutron and proton spins does not appear. The AFM 
spin order parameter arises at density g k, S.S^iq and is 
saturated at « 4.3£io- These results may be of impor- 
tance for the adequate description of spin related phe- 
nomena in the interior of neutron stars. 

A.I. is grateful for the support of the Topical Program 
of APCTP. 
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